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1 Introduction. 

Surprisingly, four dimensional gravity admits an integrable sector. It corresponds to ansatz metrics which 
admit two surface orthogonal Killing vectors. They describe stationary axis symmetric situations or colliding 
gravitational waves depending on the nature of the Killing vectors. For these ansatz, the vacuum Einstein's 
equations reduce to the Ernst's equations m. An infinite dimensional solution generating group for them 
was constructed by Geroch [Q and later identified by Julia [^ as the afiine SL{2, R) Kac-Moody group. The 
integrable character of Ernst's equations was deciphered by Belinskii-Zakharov in ref. B. See also ref.0 for 
a discussion relating these two aspects. Since then, various methods have been applied to solve 2d reduced 
Einstein's equations: either using Backlund or solution generating transformations or the Belinskii-Zakharov 
method, cf eg. g, 0, g, g, |^ [ll[ for a sample of references and ||l^, |3[ and references therein, or using 
analytical finite-gap techniques |14| . However, the by-now standard vertex operator approach to integrable 
models developed by the Kyoto school [Q was never applied to this problem. This is probably due to the fact 
that the Belinskii-Zakharov method involves the so-called moving poles which forbids a direct application 
of the vertex operator technique. This problem was partially overcome in ref . pq| . The aim of this paper is 
to fill the gap left open in this work and to describe how vertex operators may be used to solve Einstein's 
equations. It leads to determinant formula for the metrics which are described in the following section. This 
vertex operator approach, which is based on an algebraic formulation of the dressing group, may also be 
useful for deciphering the Lie-Poisson properties of the solution generating groups of Einstein's equations. 
See refs.[^ |l8) for a discussion of the Lie-Poisson properties of the dressing transformations and ref. [^ for 
a discussion of these properties for the Geroch group. In supergravity context these groups are called duality 
groups and they are important for quantization. 

Besides providing explicit formulas for exact solutions of Ernst's equations, (but whose possible physical 
applications are not discussed) , one of the aim of this paper is to decipher the algebraic structures underlying 
the solvability of Ernst's equation. Contrary to the impression that one may get from the Belinskii-Zakharov 
approach which uses space-time dependent spectral parameters, we shall show that Ernst's equations belong 
to the usual class of integrable hierarchies, such as the KP or sine-Gordon equations, and that they may be 
solved using usual algebraic technical tools. Vertex operators is one of those techniques which provide tools 
for algebraically solving matrix Riemann-Hilbert problems. Deciphering these algebraic structures leads 
us to identify dual pairs of elements of an afhne Kac-Moody group corresponding to solutions paired via 
the Kramers-Neugebauer duality. This could be interesting in view of the importance of the duality group 
in supergravity analogues of Ernst's equations. But the long term motivation for this algebraic detours is 
quantization. 2D reduced gravity provides a toy model in application to quantization technique to gravity. 
The dressing method, on which our approach is based, is particularly adapted to quantization since the 
dressing group usually acts on the phase space of integrable hicrachies by a Lie-Poison action. So it is 
promotted to a quantum group symmetry after quantization. The next step in that direction would thus be 
to find a concise description of the symplectic struture of this system and of the Lie-Poisson property of the 
dressing group. 

We have tried to write the paper such that it may be read in two different ways depending if one is only 
interested in explicit formulas for the solutions or if one is willing to learn the algebraic structures underlying 
the derivation of these formulas. Readers who are just interested in concrete formulas to obtain solutions 
may restrict their attention to Section 2, beginning of Section 3 and Section 6. Algebraically oriented readers 
may look at Section 4 and Section 5. The main new trick which allows us to complete the approach initiated 
in pa is a construction of dual pairs of elements in the SL{2, R) Kac-Moody group associated to dual pairs of 
solutions exchangeable by Kramers-Neugebauer involution. This leads us to solve a problem of factorization 
that was left opened in p6| and to give determinant formulas for the metrics. 

The vertex operator and the determinant formula of the metrics are described in the following Section 2. 
The rest of the paper is devoted to the proof of these formulas and it is organized as follows. In Section 3, 
we recall basic facts concerning the 2d reduced Einstein's equations and present them in a way convenient 
for the following. In particular we introduce the appropriate tau functions. Section 4 is devoted to the 
construction of dual pairs of vertex operators which allow us to compute the dual pairs of tau functions. 
Section 5 describe how to get the coefficients of the metric given the dual pair of vertex operators. Finally, 
a few explicit examples and comparisons with the previous result are described in Section 6. Appendix 



A presents a rapid survey of the method described in |16 , and we have gathered a few useful formulas in 
Appendices B and C. 



2 Determinant formula for the metrics. 



As is well known, solutions of the reduced Einstein's equations come in pairs which are related via the Kramer- 
Neugebauer duality. The two dual metrics, that we shall denote by ds^ and ds^, can be parametrized in 
terms of Weyl coordinates z and p as: 



ds^ 
dsl 



p2CT /j,2 
p2ct* fj„2 



p 5 e'" {dz^ - dp^) + Gabdx'^dx'' 



p-^e'" {dz^-dp^) + Gli,dx''dx'' 



(1) 
(2) 



All fields only depend on the two coordinates z, p. The indices a, b run from one to two and p^ = detGab = 
detG'*j,. The prefactor ct, or (?*, is usually called the conformal factor. The precise form of the duality 
relation mapping ds^ into dsj is recalled in the following section as well as an alternative parametrization 
of the metrics. 

The components of the metrics would be parametrized in terms of expectation values of certain vertex 
operators as follows: 



s'"Gi2 = Vp Im 



e'^'=l(5*>..pl 



e^'^Gn 






o2ct* 



„2ct* 



G: 



22 



VpM. 



G*i2 = \fP Im (g)z^p {Wg):,^ 



iz,p\ 



e'^ Gl, = ^p\{Wg),.J 



(3) 
(4) 
(5) 
(6) 



Here g, Wg and g* , Wg* denote the vertex operators associated to the two dual solutions. See eqs.( [42| , |43D 
and eqs.( |5qj59| ) below for their definitions. The overbar denotes complex conjugation. The indices z, p are 
here to recall that these expectation values depend on the Weyl coordinates. Notice the interplay between g 
and g* in the above formula: the element g enters in the conformal factor a while the dual element g* enters 
in the metric Gab- 

The solutions we shall describe depend on two sets of real parameters with total number 2{m + n): the 
first set is made of 2m parameters (zp,Up), p — 1, • • • ,771, while the second set is made of 2n parameters 
denoted {zj,yj), j — 1,- ■ ■ ,n. All parameters are real but \zj\ > 1 and \zp\ > 1. It is convenient to introduce 
the notations Xj, pj and Xp, pp such that: 



X,=p 



z]-\ 



— — a,iid //,- — 

(Zj — Zy + p^ ■' (2j - z) - p 



(7) 



and similarly for Xp and Pp. The functions pj are related in a simple way to the usual moving poles in the 
Belinskii-Zakharov approach. The validity of the metric is restricted to (z, p) domains such that z ± p < 2: j 
for Zj > 1 and z :L p> z^ for Zj < —1, and similarly for Zp. 

The explicit determinant formula for the metrics come from the following expressions for the vertex 
operator expectation values which we shall derive in the following sections: 

(8) 
,) (9) 



(10) 



{g).,p^n ■ t{Y,\p,) ; 




(g*).,p = pi r!* • r{Y,B,\p, 


{Wg)..^p^pn^ ■ TiY,B^\p,) ; 


{Wg*).,p ^p--ni ■ t{Y, Br 


The tau functions t{Y\p) could be written as n x n determinants: 


t{Y\p) = det„xn 


l + iV 




with y, - 2^^^^- 
Pi + Pj 


The parameters are the following: 




/ m \ 

Y.^yjX, 11 b;; with B,p 
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^^'^-^^ and B, = ^-^'^ 
Pj + Pp 1 + P] 


Finally the prefactors fi and 17^, are given by: 
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The coefficient Gab of the metrics are thus the ratio of these determinants. Note that the prefactor fi cancels 
when computing these ratios. 

Stationary axis symmetric solutions may be obtained by analytic continuation, see Section 6.3 below. 
However the reality conditions are then more involved and not all of these solutions are physical. Notice also 
that for an infinite number of parameters {zj,yj) the above solutions may be described in terms of Frcdholm 
determinants. 

Of course many solutions of Ernst equations have already been described in the literature, see P, R, pi 
P, 11^, 0| for a sample of references and O, 0| and references therein, The formula closest to the ones 
we are describing here are those obtained by P. Letelier, cf. ref.[|lO|. These later formulas also present the 
metric coefficients as ratio of determinants which are analogous to those we have described but are not quite 
the same. Moreover in the formula [[lO| a few quadratures, which in the formulas above are absent, remain 
to be solved. But, again our principal aim was to show explicitly how vertex operator methods, which are 
standard in integrable models, may be applied to the Ernst equations. 

The formulas ( |3|-p[ ) for the metrics originate from an algebraic expression for the Ernst potential which 
was derived in ref. |]l6| using the algebraic method of dressing transformations: 

G22 'gI, (vl/o • (g-ig+) • *-!> (g).,p ^ ^ 

Here, *o, M/2(l) and g± are operators acting on an auxihary Fock space and (■ ■ ■) denotes the vacuum 



expectation value in that Fock space. ^0 is an explicit operator defined in eq.(|32[) below which carries all 
{z,p) dependence. ^2(1) is a specific vertex operator associated to the Kramer- Neugebauer duahty and 
g = gZ^g+ are elements of an SL{2, R) affine Kac-Moody group. 

However, the formula (|l3| ) is not the final answer since it requires factorizing the element g in two pieces 
as g = gZ g+- This factorization problem is specified by Einstein's equations up to a residual 5*0(2) gauge 
freedom which refiects the gauge symmetry of these equations. Eq.(|l^) is only true in a specific gauge: the 
triangular gauge. The problem of fixing this gauge in the factorization of the group element as 5 = gZ^g+ 
remained unsolved in |16|| . The main new technical point of this paper is to solve algebraically this gauge 
condition. This will be done in three steps: (i) by relating the triangular gauge condition to the dual pair of 
solution of Ernst equations, cf Proposition 1, (ii) by establishing the relation between pairs of group elements 
g and g* corresponding to dual pairs of solutions, cf Proposition 2, and (iii) by using these two results to 
factorize the element g ~ gZ^g+ and to compute gZ^W2{i)g+, cf Proposition 3. This then leads to the 
determinant formula (§,^). The key algebraic result is the duality {g,g*), see eqs.(^,^, or eqs.(62 6^) for 
its interpretation in the dressing group. 

Although this paper is necessarily quite technical, we have tried to minimize the technical aspects. In 
particular the proofs could probably be omitted when reading it. Since this paper is a continuation of ref . [[l6[ 
and has for aim to fill the gap left open there, we shall freely use results from that reference. 

3 Equations of motion and dualisation. 

This section is devoted to recall a few well known facts on the Ernst equations: its duality property and the 
relation between this property and the triangular gauge condition. 
Following Papapetrou, let us parametrize the metric (j^ as: 

ds^ = p~5e2^ (dz^ - dp^) + pA-ida;2 _,_ ^^(^^ _ j^^^y ^^4) 

This means that G22 = pA and G12 — —pAN. The fields (a. A, A'') only depends on the two coordinates 
z and p. This ansatz can be made more covariant by introducing two light-cone coordinates u and v such 
that 

p = a{u) + b{v) , z — a{u) — b{v) (15) 

with a and b any functions. The metric is then: 

ds^ = -Ap-ie^" dadb + pA'^dx^ + pA{dy - Ndxf 

Consistency of this parametrization follows from the fact that Einstein equations imply that p defined as the 
square root of detGab is harmonic: dudvp = 0. In the following we shall denote by d+ and d- the derivatives 
with respect to u and v. 



As is well known, given any solution of the reduced Einstein's equations one gets using the Kramer- 
Neugebauer symmetry a dual solution by introducing the metric ds^ parametrized as above but with the 
original fields (ct, A, N) replaced by the dual fields {a*, A*,N*) with: 

A* = l/(pA) , A*d^N* = e^aA^aiV , A*e^^' = Ae*^ (f6) 

Hence, the dual metric dsl can be written as: 

dsl ^ Ae^'^ {dz^ - dp^) + p^Adx^ + A-^{dy - N*dxf (17) 

3.1 Reduced vacuum Einstein's equations. 

The vacuum Einstein's equations, which codes for the Ricci flatness of the metric (nh, reduce to the Ernst 
equations. They are those of a non linear sigma model on the coset space SL{2,R)/SO{2). They can be 
written using a first order formalism. This amounts to introduce a connection taking values in the sl{2,R) 
Lie algebra but covariant with respect to the so(2) subalgebra. We denote by Q and P its components with 
Q antisymmetric and P traceless symmetric 2x2 matrices. The equations for Q and P are then, cf eg. 
ref. ^21 s-nd references therein: 

[d+ + Q+ + P+,d- + Q- + P-] = (18) 

D-{pP+)+D+{pP-) = (19) 

2id±d){p-'d±p) - tr{P±P±) (20) 

where D± = d± + [Q±, ■ ] is the covariant derivative. This system is gauge covariant with gauge group 
SO{2). It acts as Q± ~^ Q±+ A~i9±A and P± —f A~^P±A for A G SO{2) while a is gauge invariant. 

The first equation ( [l8| ) means that {Q± + P±) is flat; this ensures that there exists a 2 x 2 matrix V, 
with V € SL{2, R) such that Vd±V^^ = Q± + P±. The matrix V is the matrix of the zwei-beins such that 
the metric coefficients Gab are such that: 

Gab = p{'VV)^^ 
Note that detGafc = p^ since detV =1. It is often convenient to choose a gauge in which V is triangular: 

/ A-i \ [A-i \ / 1 0~ 
\-NAi Ai J \ Ai J \-N I 

It reproduces the parametrization (04) for the metric. In this gauge the connection is: 

Q±+P± = i(A-i9±A) • a' + A{d±N) ■ a~ (21) 

Plugging this parametrization into cqs.(18 1^) gives the Ernst equations. 

3.2 Alternative parametrization of Ernst equations. 

To make contact with the vertex operator approach described in the following we need to introduce an 
alternative parametrization of the connection. Since Q is antisymmetric and P traceless symmetric, we may 
parametrize them as follows: 

Q± ^ I {d±cf>±) [a+ - a-] (22) 

P± = ^ip-'d±p){Z^e'^^) [a'-i{cT++a-)]+c.c. 

with CT^, a Pauli matrices. The dynamical fields are thus the two real fields (/)± and the two complex ones 
Z±. Gauge transformations act by translations as (j)± — > (l)± + X with A a real function. In particular Z± as 
well as = 0+ — (f)- are gauge invariant. 



In this parametrization, the Ernst equations become: 



d±Z± = \{p-^d±p){Z±-e^'^Z^) 

d+d-cf - -'-{p-'d+p)ip-'d-p) {Z+Z^e"^-Z^Z_e- 
45±a = ip-'d±p)Z^'Z^ 



(23) 

(24) 
(25) 



These equations are gauge invariant since (f> and Z± are gauge invariant. The equivalence between eqs.(p3[) 
and the more usual Ernst equations follows by plugging the parametrization ( |22| ) in the first order formalism 
eqs.(|8|l9). Note the similarity with the sine-Gordon equation. 

Let us also mention that eqs.( p3|p4|j25| ) may be rewritten in a bilinear form similar to Hirota's equations. 
Namely, let tq and t± be tau-functions defined by: 



To = exp(cr - -c, 



T+ = Z+ To, T_ = Z_ To 



then, eqs. (|23| , p4| , p5|) are equivalent to the following bilinear equations: 

1 



■^{P ^5±p) (ror± - roTzp) 



ro (9±T±) - T± {d±To) 

To (9+5- To) - (5- To) (5+ To) 

This is proved by direct substitution. Eqs.(E6h are Hirota's form of Ernst's equations. 



'^{P ^d+p){p ^d-p)T+T- 



(26) 



3.3 Dualisation and triangular gauge. 

Let us now explain the interplay between the condition for choosing the triangular gauge and the existence 
of two dual solutions. 

Proposition 1. Let {(t)±,Z±) and {(p'^jZ'^) be the components of the two dual solutions assuming that 

both connections are in the triangular gauge. The duality relation is then equivalent to: 

— a duality relation between the helds (f>± and 0j_ , which is valid only in the triangular gauge: 



'■^-\ 



±i 



±Q 



'±' 



a duality relation between Z± and ZZ^ , which is gauge invariant: 

ZI + Z+ = -e-'(0+^*)/2 
Z_ 



(27) 



(28) 



YZ = -e-^'^'^-^")!'^ 



Proof . Indeed, let ((/)±, Z±) be a solution of the Ernst equations (23,2J). Imposing the triangular gauge as in 
eq.(pi|) demands that the coefficient of Q± + -P± along a^ be zero. In the parametrization (^2|) this amounts 
to impose: 

2a±</.± =z(p~i9±p) {Z^e'^^ -^e-*^T) (29) 

Let us now define the dual solution (/>* and Z|. by the formulas (p^) and (ESh . Notice that when parametrized 
in terms of and 0*, the condition (p9|) becomes: 



a± (0* ± (/.) = i{p-^d±p) 



Z^e^ 



(0*T0) 



Zz^e^s 



('t>''T4>) 



One has to prove that (0*, Z^) are solutions of eqs.(23,2J). Consider eq.(p3|) for Z^. Using the definition of 
Z\_ and eq.(|^ satisfied by Z+ and eq.(29) for + 0*, we get: 



d+Zl 



1 



{p-^d+p) 



\{p-'d+p) 



Z_ 



Z*_e-''> 



Z_e2 



Z_p/ 



e 2 



(0'+0) 



This shows that Z^ is the solution. The equation for Z*_ is proved in a similar way. The equation for the 
dual field 0* is proved by taking derivatives of cq.d29). Finally, since the transformation (pi) is an involution 
the triangular gauge condition for 0j_ follows from that of 0, eq. (p9|) . ■ 



Remark that the duality relation (G^) and (pa) are purely algebraic and that eq.(E7h is similar to the 
usual abelian T-duality of string theory. This proposition also shows that solving for the triangular gauge 
or for the dual solution are equivalent problems. 

Using the defining relation between 0, Z± and the connection Q± + P± in the triangular gauge, one may 
translate the duality (p^pSl) on the connection as: 

Q± + P± = i(A-i9±A)-a^+A(a±iV)-a- (30) 

Q\ + Pl = -i(p-ia±p)-a--i(A-i9±A).a^±A(9±7V).a- 
This is clearly equivalent to the Kramer-Neugebauer duality (|lq). 

4 Integrability and vertex operators. 

The aim of this section is to recall the use of vertex operators for solving the Ernst equations following the 
method initiated in ref.[|l^. This method was based on an application of the dressing method [|^. We 
have written this section in a down-to-earth way by just presenting propositions which give the rules for 
computing solutions of the Ernst equation using vertex operators. These rules are analogous of the famous 
Kyoto formula flq] for the tau function of the KP hierarchy in terms of fermions or vertex operators. The 
logic which allows us to go from the dressing method to these vertex operator formulas is recalled in Appendix 
A. We shall also describe the dual pairs of vertex operators corresponding to dual pairs of solutions 

4.1 Vertex operators and tau functions. 

Vertex operators are exponentials of a free bosonic field acting on an auxiliary Fock space. They may be 
used to find solutions of the Ernst equations in a way similar to their use for solving the KP hierarchy as 
described the Kyoto school flS]. Let us denote by Ar(w) the bosonic field: 



X{w) = -i V" p_„ with [p„,Pm] = n5„+m,o 



n Oi 



dd 



The operators p„ generate a Fock space, we denote by |0) its vacuum: p„|0) = for n > 0. For any number 
u, let Wu{w) be the vertex operators: 

Wuiw) =: exp{-iuX{w)) : (31) 

The double dots refer to the normal ordering which amounts to move to the right the oscillators p„ with n 
positive. The parameter w is called the spectral parameter. The Virasoro algebra acts on the Fock space 
generated by the j)„. The Virasoro generators L„ are represented by 



Y^iLn ~ ^A..o)w-'"-' - - J : {du.Xf : 



16 ' ' 4 

n 

The algebraic dressing method applied to the Ernst equations leads to the following result: 
Proposition |l6| . Let ^o be defined as: 

Let g he any product of vertex operators of the following form: 

m n 

5 = const. X{Wu^(wp)-\{{l + iyjW2{wj)) (33) 



where {yj, Wj) and {up, Wp) are 2(rn + n) real parameters, then the fields [a, 0, Z±), defined by the foUowing 
expectation values, 

exp(^a-^0^ = (*o5*(7') = (.g).,p (34) 

Z7 exp (d ~ ^0 j = ((1'o.g*(7^) • P-i) = {9P-i)z,p (35) 



Z_ exp ( a - -0 ) = -(pi • (*o5^o"')) = -{pig)z,p (36) 



4' 
are solutions of equations ( [2^j24|^ - 

The above equations serve as the definition of {g)z,p, {gP-i)z,p and {pig)z,p- This is the notation used in 
Section 2, eqs.(|-g|). These expectation values are the tau-functions of the model. 

To make sense the vertex operators in eq.(|3^) have to be ordered in decreasing order of the Iwjl's. 
Modification of this order may be done by analytic continuation in the spectral parameters. 

To compute these expectation values one needs to know how to conjugate vertex operators Wu{w) with 
*o- One has @: 

*o • Wu{wj) ■ *o ' = xf'^ ■ Wu (m,) • (37) 

with Xj and /ij defined in eq.(0) with the parameter Zj = Ji^i - Notice in particular that ^Pq ■ W^u(l) ■ ^o ^ = 



/o^r Wu{^)- With this result in hand, the computation of the expectation values ( p4p5y3q ) is reduced to the 
computation of expectation values of vertex operators. As recalled in Appendix B, this is done using the 
usual Wick's theorem. For example: 

ra ni n 

{g)z,p = (11 ^p'") (IT ^-M) • n (1 + *^.^2(/i,))) 

p=l p=l j=l 

= ^ ■ r{Y,\^l,) (38) 

with fi and T{Y\y^) defined in eqs.( pi]| , p^ ) above. The expectation values {gp-i)z,p and {pig)z,p are computed 
similarly using the formula recalled in Appendix B. They may be obtained from the previous expression for 
{g)z,p by replacing each monome Xk^ ■ ■ ■ Xk^ by, 

Xk,---Xk^ -^ Xk,---Xk^ i^Uk,^^J■k^] for {gp-i)z..p (39) 

^fei • • • ^fcp ^ - Xk^--- Xkj, I ^ Uk^fik^ I for {pig)z,p 



in the formula (38) 



4.2 Vertex operators and duality. 

Two dual solutions correspond to two dual products of vertex operators that we shall denote by g and g* 



The duality relations (E7p and (28) can be translated into quadratic relations for the expectation values of 



these operators similar to the Hirota equations. Namely: 



{g)z,p {g*p-i)z,p + {g*)z,p {gp-i)z,p = -{g)z.p {g*)z,p (40) 

{g)z,p {pig*)z,p + {g*)z,p {pig)z,p = +{g)z.,p {9*)z,p (41) 

Proposition 2. Pairs of solutions of the above duality relations, eqs.( [ic| , [il| ), are provided by the following 
pairs of vertex operators: 

m n 

g = const. l[Wu^{wp)-l[{l + iyjW2{wj)) (42) 

m n 

g* =. const. W-iil) ■ J] W^u,{wp) • J] (^ + WjW-2{wj)) (43) 

p=i j=i 

8 



Up the multiplication by W-i{l) the dual vertex operator is obtained by the charge u into — u. The constant 
pref actors in eqs. m2S43) are irrelevant. 



Proof . It relies on an identity for the tau- functions proved in ref. p2] . Let us sketch the proof of eq. ( fJOJ ) for 
m = 0. Recall eq.(|8|) for {g)z,p and eq.(^) below for {g*)z,p- 

{9)z.p--T(Yj\fi,) and (g*)z,p ^ p^ T{BjYj\fij) (44) 

From eqs.py) one infers that the expectation values {gp-i)z,p and {g*P-i)z,p rnay be written in terms of 
derivatives of tau- functions. Namely: 

d 

{9P-i)z,P = —T{Bojyj\Pj)\po=o 

P'^ {g*P-i)z,p = -t{Bj Yj\pj) - —TiBojBjYj\pj)\f,„^o 



by adding and expending in power of /iq the following two relations proved in 



with Baj ~ _ ^ . Eq.(|40|) may then be written as a bilinear identity for the tau- functions. The later follows 



T{BokBikYk)T{Yk) + T{BokBikYk)T{Yk) = T{B^kYk)T{BikYk) + T{BokYk)T{BikYk) (45) 
Sio [T{BakBikYk)T^) ~ T{BokBikYk)T{Yk)\ = r{BnkYk)T{BikYk) - T{BokYk)T{BikYk) (46) 

with Bik = -Bfc and Bio = \+^° ■ "^^^ proofs of the general case m 7^ as well as eq.(^) are similar. ■ 

In eq.([43|) the order of the operators matters: the operator VF_i(l) has to be on the left. When changing 
the order of the operators one has to take into account their commutation relations, ie. IV_i(l) anticommutes 
with W-2{w). 

Expectation values of the dual vertex operators may be evaluated using the conjugation formula (|37 
and Wick's theroem as explained for {g)z,p- One gets: 

{g*)z,p ^ {^09*^0') 

- pin* ■ TiB,Y,\^i,) (47) 

as in eq.(||). 

5 Algebraic computation of the metric coefficients. 

The previous section made precise the relation between the two dual vertex operators g and g*. It allows 
us to compute algebraically the gauge invariant fields (a,(j),Z±) and their duals {a* , cf)* , Z^.). This is not 
quite the final answer for the metric since to go from {(f), Z±) to the metric coefficients Gat, or A and TV, one 
needs to impose the triangular gauge and then integrate the connection Q± + P± to obtain the zwei-bein V. 
Imposing the triangular gauge requires solving eq. ( p9| ) . This is a non- linear problem which was actually not 
solved in [|l6|. We shall now solve it using our knowledge on the dual pairs of vertex operators, eqs.( [4^ , [43|) , 
and on the link between the duality and the triangular gauge, eqs.(^,^8|). We will then be able to use 
formula ( ^3] ) to compute the metric coefficients Gab- 

5.1 Vertex operator representation and factorization. 

In order to be able to apply formula dT3) we need to make a small detour into group theory in order to 
explain the relation between factorization problem in affine Kac-Moody group and vertex operators. 

This relation arises because the vertex operators ( pl| ) may be used to represent the sl{2, R) affine algebra 
on the Fock space. The commutation relations of the sl{2, R) affine Kac-Moody algebra are: 

[x®r,y®n = [x,y]<E,t^+"' + n^tr{xy)5n+m,o (48) 

The affine Kac-Moody sl{2, R) algebra is twisted in the sense that its elements x^t" are such that x E so{2) 
if n is even while x is an 2 x 2 traceless symmetric matrix if n is odd. We are actually considering the semi 



direct product of the Virasoro algebra with the affine algebra. The crossed Lie bracket is: [Ln,x(E)t'^ 

The representation of the affine algebra on the Fock space is specified by the following relations [ES^ 

. dXiw) 



dw 



: Y^{cj'-® r)w-" (49) 

n odd 



y ((a+-a-)0r)^-"-2y 

n even n odd 



This is a highest weight representation. The highest weight vector is identified with the vacuum vector |0). 
It is such that [a^ - ct~)|0) = ||0) and {(t+ - a~) (g) i"|0) = for n > 0. Remark that in particular eq.(|9|) 
means that (a^ (g) i") is represented by the bosonic oscillator p„. 

The algebraic dressing method from which equation (|l3| ) follows relies on a factorization problem specified 
on the affine SL{2, R) Kac-Moody group. This is defined as follows. Let B± be the two Borel subalgebras 
respectively generated by the central charge k and the elements x (g i^" with n positive. Let B± — exp B± 
the corresponding Borel subgroups. Then the factorization amounts to formally decompose any element g 
in the affine SL{2, R) Kac-Moody group as the product of elements in B±: 
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= gjg+ with g±eB± (50) 



where we also demand that the components oi g± on the exponential of the central charge is inverse. Remark 
that since the two Borel subalgebra B± have elements x (g i" in common, the above factorization is defined 
only up to a multiplication g± — > hg± by elements h G SO{2). As explained in ref . |Lq| , this freedom is linked 
to the SO (2) gauge symmetry of Ernst equations ( [lq , [l^ ). 
We can now state a result from ||T^ : 

Proposition [p^. Let £ = ^+(1) — ^-(1), let g — gZ^g+ be factorized according to the triangular gauge, 
then 

A* (vl/o • (ffZ V) • *o-') 

Furthermore, in the Fock space the elements £ is represented: £ = iW2{l)- 
A similar expression holds with the dual operator g* and A and N . 

Eq.(|5l|) is only valid in the triangular gauge. Let us explain in more detail how this gauge choice fixes 
the factorization of g as gZ g+. Since g± G B±, one has: 

g± = exp(±-kj exp(^-^(cr+-cr")j x [| degree] > 1] 

For a given group element g = gZ g+ only the difference ip^ — ip^ is fixed. To translate ip± keeping this 
difference fixed amounts to multiply g± —^ hg± with h S S0{2). The link between this freedom and the 
SO{2) gauge freedom of the Ernst equation relies on the fact iQ that ip± coincide with the fields (j)± at 
z = Q, p = 1. Thus in the triangular gauge, 

V±^\{V*±V) (52) 

with if and (p* equal to <& and <&* at z = 0, p = 1. This follows from eq. (pTj) . Since ^0 = 1 at z = 0, p = 1, 
they may be evaluated using eq. (p4|) : 

-(■1)^1 • -(■f)^i <-' 

Here, the expectation values are the vacuum expectation values, without insertions of $o- 
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5.2 Factorization and dualisation. 

The factorization problem in the afhne Kac-Moody group may be understood as a kind of normal ordering. 
So when considering the vertex operator representation one has to face two different normal orderings: the 
one associated to bosonic oscillators p„ and the group theoretical one. We shall now explain the link between 
them for the vertex operators ([4^). 

First, consider vertex operators Wu{w). Since {a^ (^ f") is represented by p„, we may consider them as 
elements of the affine group. Namely: 

W^^H = -p(-^.^«log(i±^)).exp(|.^«log(i±|^)) (54) 

= W^{w)z'-Wu{w)+ 

The last equation serves as definition oiWu{w)±, elements of the Borel subgroups B±. Thus, the two normal 
orderings coincide for these group elements. 

Consider now the product of vertex operators of the form J^ ■ (1 + iyjW2{wj)). Since W2{'w), which are 
generating functions representing elements of the affine algebra, are nilpotent inside any correlation functions, 
ie. W2{w)W2{w) = 0, these products are representations of elements of the Kac-Moody group. As shown in 
ref. |2^ these products may be factorized in the affine Kac-Moody group. More precisely, let g±{j) be the 
elements of the Borel subgroups B± defined by 

9±{j) = exp (±^fc) exp (^(a+ - a')) exp (y4(u;,)) (55) 



with 



£±{w)^± 



l + (i/it;)±2\ _ / 2(t/«;)±i 



(^ - '^ ) ® 1 //...+, ]-i^^+ ^') 



Then, for fc = 1, • • • , n one has 1221, 



(56) 



gZ\l) ■ ■ ■ gZ\k) ■ g+ik) ■ • •5+(l) = l[{l + iyjW2iw,)) (57) 

Eq.(p^) is valid in the Fock space representation. The relation between the parameters {sj, rj, Vj) and {yj, Wj) 
is explained in the following proposition. 

We can then solve for the factorization problem: 



Proposition 3. For g and g* the dual vertex operators (63) and {43), then: 

m n 

Wg = gZ'W2{l)g+= {-la + bW2{l)) -Jl Wu^{wp) -Hil + iyjW2{wj)) (58) 

7n n 

Wg* = g*_-'W2{l)gl = {-la* +b*Wi{l)) -Hw^uAwp) -Hil + iy,W-2{w,)) (59) 

in the triangular gauge. Since the Ernst potential is defined up to a multiplicative constant and a constant 
translation on N, a and b are irrelevant when computing the metrics. 

Proof . This relies on the relation |22J between the parameters {yj , Wj ) involved in the vertex operators and 
the parameters isj,rj,Vj) in the group elements g±{j) such that: 

n 

gZHl) ■ ■■gZ\n) ■ g+{n) ■ ■ ■ g+{l) - \{ {\ + iy,W2{wj)) (60) 

Let us introduce the obvious notation Lpj and Lp* for j < nhy 

j j 

fc=l k=l 
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such that Sj = —\{ipj — (fij-i) and similarly for Vj. Of course ip± = ^{(p^ ± 'fin)- First {sj,rj) are recursively 
determined as functions of {yk,Wk) with A: < j by computing the expectation values of eq.(p7|): 

exp(^(rfe+z|i)] ={f{{l + iykW2{wk))) 

\k=l I k=\ 

These equations determine fj{yk) as functions oi yk, k < j. The fj's are then given by |2^ : 

Wj+l - Wk 



y^jiVk) = ^jif3j+i,kyk) with /3j+i,fc 



Wj+i + Wk 



This leaves Vn, which actually cancels in eq.(60), undertermined. However, w„, or equivalently (^*, is fixed 
once we impose the triangular gauge. Indeed the triangular gauge condition written as in eq.(p3|) and 
formulas (^ 43) for the dual vertex operators leads to: 



VniVs) = Vnitijyi) with Pj = — — i (61) 

i -f- Wj 

This allows us to go one step further in the recursion relation ( |57| ) by inserting one extra vertex operator 
with spectral parameter w equals to 1. Let g±{n + 1) be such that gZ^{n + l)g+{n + 1) = a + ibW2{l) with 
a, b numbers. The triangular gauge condition (pl|) then implies: 

n 

9Z\l) ■ ■ ■ gZ\n) ■ {a + ibW^il)) ■ g+{n) ■ ■ ■ g+{l) = (a + zbW2il)) l[{l + WjW^iw,)) 

with a,b functions of a, 6. Taking a — and 6=1, this proves eq.(p8|). The dual equation ( |59| ) is proved 
similarly. ■ 

Once the operators gZ^W2{l)g+ and gZ^^W2{l)g'*l are expressed in terms of the product of vertex 
operators, it easy to conjugate them with ^I'o using eq.pTI). One may then evaluate {Wg)z,p and {Wg*)z,p 
using Wick's theorem, cf. Appendix B. Of course one gets eq.(||). This ends the algebraic proof of the 
determinant formula for the metrics. 

Remark that eq.([45|) for the tau- functions with /^o = 1 implies: 



Re((g).,p(W^5).,p)=N/p|(5* 



'z,pl^ 



This shows that G*^ defined in cqs.(0-@) satisfies detG*^ — P^ ■ It provides a non-trivial check of the 
construction. 

5.3 Dualisation in the dressing group. 

The dressing group is the group whose elements are the pairs {g-,g+) factorizing g a.s g = gZ g+- It is 
different from the afhne SL{2, R) Kac-Moody group since their multiplication laws do not coincide, |1^. In 
the dressing group the product is given by {g^,g+){h^, h^) = {g^h^,g+h+). 

The solutions we have obtained should actually be labeled by elements of the dressing group since this 
is the solution generating group. As a consequence, the duality between the vertex operators (E3) and MS) 
should be thought of as a duality in the dressing group. Writing the vertex operator (B2|) in the dressing 
group amounts factorizing them according to the rules explained in the previous section: 

m rn 

9-^9+ = n WuA^p)Z' ■ gZ\l) ■ ■■gZ\n) ■ g+{n) ■ ■ -g+il) ■ J] Wu,{wp)+ (62) 

p— 1 p— 1 

where the middle term corresponds the factorization of Y[j=i (1 + ^Vj ^2{'Wj)) according to eq.(p7[), 

n 
gZ\l) ■ ■■gZ\n) ■ g+{n) ■ • •g+(l) = Y[{l + ty, W2{w,)) 
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Here yj — yj Y\ Ap ■ Similarly the dual vertex operator is factorized as: 

m 



p=l 



where 



ff- 



•?;(n) •• •?;(!) • n W^-«.K)+ ■ W-iil)+ (63) 



-i(i)...5r'N-?;(n)...g;(i) = n(i + *^;^-2K-)) 



with y* = ft^?^ 

In eqs .( 1621 , 1631 ) the dual elements {g^,g+) and {g*^,gX) are written as elements of the dressing group (and 
not in a particular representation). It is then clear that the map from {g-,g+) to {g*_,g'^) is an involution. 
However it is not a group automorphism R It is unfortunate and frustrating that we do not know how to 
write this involution in more group theoretical way without relying on these particular elements. A better 
group theoretical understanding of the dualisation, ie. of the relation between g and g* , will provide a way 
to decipher how general the duality property and the involution trick are and whether they apply to other 
integrable systems. 

6 A few examples. 

Since our aim was to describe the use of vertex operators for solving Ernst equations and not to describe the 
physical properties of the solutions, we shall only discuss a few examples (which are actually already known 
in the literature, cf eg. |^, 0, g, g| and Q and references therein). 

6.1 Diagonal solutions. 

Diagonal solutions correspond to solutions for which Gab is diagonal, ie. A^ = 0. They are obtained by 
imposing yj = in the parametrization of Section 2. They correspond to vertex operators of the form 

ni 

gZ'g+^Ww^^iwp) (64) 

p=i 

They depend on the 2m parameters (zp, Up); recal that Zp — l_. . Here and in the following examples, we 
drop insignificant multiplicative constants in front of the vertex operators. 

The simplest of such solutions is the well known Khan-Penrose metric Q . It describes the interaction 
region of two plane impulsive gravitational waves having their polarization vectors aligned. For more details 
on this subject, see eg ref.[^. It corresponds to the following set of parameters, {(zp, Up)} = {(1, 1), (—1, 1)}, 
or equivalently to the following vertex operator: 

Wi{oo)Wi{0) 

Notice that this choice of singular values for Zp leads to null values for Xp. However, we may reabsorb these 
singular constant factors in the normalisation since Einstein's equations determine the conformal factor only 
up to an additive constant. In order to have the formula closed to those which may be found in the literature, 
we introduce two new positive variables u and v defined by: 

p = 1 — u ~ V and z — v — u 



^Note however that the relation between g* and g may be written as g* = W—i{l)'T{g) with T the automorphism of sl{2, R) 
fixing so(2) and multiplying 2x2 traceless symmetric matrices by minus one. 
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The domain of validity of the metric is m > 0, w > and u'^ + v"^ < 1. Simple computations using eqs.(^,^,||) 
lead to the following line clement: 



1 + u\/\ — w^ + uv 1 — U^ , o 1 — u\/\ — v^ — v\/l 



u 



2 



+ (1 - ^2 - «2) ^ ^=^ dx2 + -== -== dy^ 



, 1 — u\/l — v'^ — uVl — u^ \ + u\/l^^v^ -\- vyl — v? , 

This is the Khan-Penrose solution. 

Another class of well-known solutions are provided by Kasner's solutions. They correspond to the pa- 
rameters {(zp,Mp)} = {(z,u)} in the limit z — > -f-cx). Equivalently they correspond to the vertex operator 

Wu{w) with ui ^ 1. With t = p~^~ and after a correct redefinition of the variables in order to absorb null 
factors and non significative constants we deduce that: 

ds^ = -dt^ + fP' dz^ + i2p2 ^2,2 j^ ^2p3 ^y2 
U — 1 1 — U l-hli 

with pi = p2 = 2 and P3 = 2 such that pj + pj + pi ^ 1 and pi + P2 + P3 = 1- 

u^ -I- 3 u^ -t- 3 u^ -|- 3 

6.2 Non-diagonal solutions. 

In this subsection, we will focus on two examples of solutions describing collisions of two impulsive gravi- 
tational waves with non-colinear polarization vectors: the Chandrasekhar-Xantopoulos solution Ig] and its 
dual, the Nuktu-Halil solution , and the more general Ernst family of solutions Q . To obtain more familiar 
and more compact results, we shall once more introduce two new variables ^ and w. They correspond to 
the prolate spheroidal coordinates. Details on the relation between these various coordinates and related 
identities may be found in Appendix C. 

The Chandrasekhar-Xantopoulos solution is generated by the following sets of parameters: 

{{zp, Up)} = {(1, -1), (+CX), -1), (-1, -1)} and {(z„ %)} = | (^1, -^^) > (-1- ^^^p)) 

where p"^ + q'^ = 1. Here p, q are simply used to parametrize j/i = —y-i- The corresponding vertex operator 
is: 



g = W^_i(oo)W^_i(l)W^_i(0)(l-z^^-^W^2(oo))(l + z^^-^M^2(0)) 

Einstein's equations allow to adjust metrics elements. In particular, wc can add a constant to the imaginary 
part of the Ernst potential or multiply it by a global constant factor. Using this freedom to translate the 
imaginary part by —jr-, we obtain for the Ernst potential: 

Lf22 (-f22 [9 )z,p 

p ((1 -K)' + cPw^) - n (-2e(l - ^') + Me - u;2) -f ^^{p\v? - e) + (1 - ^^))) 

1 — p2^2 — q'^w^ 

The conformal factor is determined by computing the expectation value {g)z.p- With A" = (1 — p^)^ + q^w^ 
and Y = 1 — p^^ — q'^w'^, the line element can be written as: 



2 Vi-e' 1-wy x\ p{p + i) y px \ p{p + i) 



This is the Chandrasekhar-Xantopoulos solution [pj written in the same form as in [n2| 
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To illustrate the duality formulas, we compute the dual of this metric. It is the Nuktu-Halil solution 
The dual vertex operator is : 



The dual Ernst potential is then given by: 



4(p+l) 

1 + pC + 'i-q'u' 



G*22 G*22 {g)z,p 1-p^-iqw 



(67) 



from which the metric can be computed. 

These solutions may be seen as members of a larger class of solutions, the so-called Ernst family of 
solutions M. These correspond to the following sets of parameters: 



{(l,-l),(+(»,-n),(-l,-l)} and 



1, 



q-q' \ (_. q + q' 

4(P + P')/'V '4(P + P') 



with p^ + q^ = 1 and p'^ + q'^ = 1. Introducing p, q and p', q' is just a convenient way of parametrizing yi 
and y_i. The vertex operators are: 
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T^_„(oo)iy_i(l)iy_i(0)(l + z^^-^I^2((^)) ( 



l + i .'^ \. W2{0) 



4(p + p 
The expression of the metric is a little lengthy. For the conformal factor we have 

e^'^ {dp' ^ dz^) = \{g),j'idp'^dz') 



1 nl 



de 



dw' 



(1 - e){ip+p'n]-^^r + b' - pfi\^p 
+ (1 - wniq + q'n^-^r + iq' qnl^^r) + 2(9' - q^ne w^) 



1-^2 1 _ ^,2 

1 — W 

1 + w 



'1 — w 



For the Ernst potential we have 



p ^ Gu ^ jWg*) 

G22 G22 {g*)z,p 



z,p 



'^ B 



with 



A 



(1-a^ 



{p+p') 



1-^ 



and 



B 



+i(l-it;^)^ 



(1-a^ 



+i{l-w^) 



{q + q)[^ 



1 — w 





)* 


{q + q') T-— 

\1 + w 



+ [p' - p) 



+ {q'-q) 



+ {P - P) 



l + w 
1 — w 



+ {q'-q) 



l + n ^ 




1-ej 




.^p+^^ 


n-l-\ 
2 


'\l-w) 





These are the solutions found in P]. The dual metrics can be obtained by changing n into 1 — n and 
exchanging (p, q) with (p', q'). This may be checked by comparing the dual vertex operators. 



6.3 Analytic continuation. 

Stationary axis symmetric solutions of the vacuum Einstein equations may formally be obtained by analytic 
continuation: 



ip , X —* iLp ^ z -^ z ^ y —> IT 
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(68) 



The (dual) metric then reads: 

dsl = -A-^{dT + ojdiff + Ae2°'(dz2 + dp^) + p^A^ (gg) 

with Lo = N* . However the reahty conditions for axis symmetric solutions are more involved. 

By analytic continuation, the Khan-Penrose solution is mapped into the Schwarzchild solution. Let's see 
what happen in the case of the Chandrasekhar-Xantopoulos solution. Using the standard parametrisation 
recalled in Appendix C, cf eg.||l^, one obtains the line element of the Kerr solution: 

2M^ds^ = - 1 5- dr^ H 5- sin^ Odrdcj) 

\ R^ J R^ 

r^+a"- ^^) sin^ 0d<p' + R' (^dr^ + d9^ 

with R^ = r'^ + a^ cos^ 9 and D — r^ — 2Mr + a^. Here the domain of validity of the metric is such that 
— (M^ — a^) < D < 0, which corresponds to the region inside the ergo-sphere. However, using the analytic 
continuation described above, eg. ([6^), the variables ^ and w become 



W = - (V(l - Z)^ + p2 - V(l + Z)2 + p2^ 



The domain is now such that D > 0. The solution thus describes the asymptotically flat exterior Kerr 
solution, which is stationary axis symmetric. 

Other physically realistic axis symmetric solutions would probably require infinite sets of parameters 
(zpyUp) and (zj^yj), ie. infinite products of vertex operators. In such a case the metric coefficients could be 
expressed in terms of Fredholm determinants. However the analysis of such cases is beyond the scope of this 
paper. 

6.4 Belinskii-Zakharov approach : one soliton case. 

To make contact with the Belinskii-Zakharov approach, we describe in more detail the vertex operator 
construction of the one soliton solution found in Q . The Belinskii-Zakharov approach |4| starts from Kasner's 
metric as a seed solution. Since these Kasner's solutions are obtained from the vertex operator with parameter 
{(zp. Up)} = {(-l-oo, u)}, we insert vertex operators whose sets of parameters are now 

{{zp,Up)} ^{{+00, u), {-I, -I)} and {{zj,yj)} ^{{-l,y)} 

They correspond to 

g = Wu{l)W-i{0) {l + iyW2m 

Given these values of the parameters Zj one has to compute the values of /ij using eq.(ra). To compare with 
ref. 01 let us introduce the same notation as this reference: 



1 + Ai-i P \ P'' 

Using this parametrization and appropriate normalization, we find using formulas of Section 1 the metric: 

ds^ = P"']---H<ir + C) ^, ^ ^ cosh((l + g)r + C)^,^,,^^, 

V(2 + l)2+p2 ^ '^ > cosh(gr + C) 

2sinh(5) cosh ((i — 0)1 — C) , ^„ , 

, , ^^'. pdxdy + ^V, ^' ^, ' p^-^^dy^ 

cosh{qr + C) cosh{qr + C) 

with e^ = Ay (choosing y positive) and u = 2q. This coincides with the one soliton solution of Belinskii and 
Zakharov. 
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So the formulas we have found can be used in the same way as those found following the Belinskii- 
Zakharov method. One may first consider a seed solution and then insert more vertex operators to generate 
multiple solitons solutions. As pointed out in Section 2, this leads to solutions similar to those found by P. 
Letelier pO| but expressed and parametrized in different ways and with no quadrature left. 

The parameters fij defined in eq. (M) correspond to the moving poles and the parameter y to the integration 
constants of Belinskii and Zakharov. The relation between the parameters of the vertex operators and the 
moving poles can be made more explicit. Namely, let Zj the positions of the vertex operators and /ij defined 
as in eq.(0) by /i| = 7-1-'' ^ then the functions 



are the moving poles of the Belinskii-Zakharov method. However the main difference between the vertex 
operator and the Belinskii-Zakharov approaches is the fact that no quadrature is needed to obtain the metric 
in the vertex operator approach. 



Acknov^rledgement: It is pleasure to thank Vincent Pasquier for useful discussions. 



7 Appendix A: A survey of the general method. 

To make the paper more self-contained, we recall here some of the results obtained in fll] . They could be 
helpful to understand the general method based on vertex operators. There are a few steps which we now 
explain. We will use notations and formulas introduced in Section 3 and Section 5. In particular we need to 
use the Kac-Moody and Virasoro algebras defined in Section 5.1, eq.(|4q) and below. 
Let us first introduce the Lax connection: 

A± = ±d±(io - i±i) +Q± + P±® t^^ T {d±d)- 

It is such that its zero curvature condition [d+ + yl+,9_ + A-] = is equivalent to d± = p~^d±p with 
d+d-p = and the reduced Einstein equations. As usual the zero curvature condition is the compatibility 
condition for an auxiliary linear system: 

(a±+A±)1' = (71) 

The solution ^ of that system is called the wave function. 

The simplest solution to Einstein's equations corresponds to Q± = P± = ct = 0. We call it the vacuum 
solution. It is easy to realise that it is associated to Minkowski's flat solution (using the dual metric and 
analytic continuation). For that solution the Lax connection is simply A± — ±d±{LQ — L±i). Its wave 
function is: 



*o(u,v) 



C2 



C4[Ci,C2) 



a{u) + 03(01,02) J \a{u) + 03(01,02) 



with p = a(u) + b(v) and z — a(u) — b(v) as in eq.(|l^) and where ci, C2 are constants depending on the initial 
conditions. In this paper, we take ci = C3 = ^ and C2 = C4 = 1 (choosing other values amounts rescaling 
and translating z and p). 

The algebraic vertex operator method is based on manipulating the wave function to generate new 
solutions from old ones. These manipulations are dressing transformations [[l7| . They were applied to 2D 
reduced Einstein equations in ref . [|l6[ . Dressing symmetries are associated to the factorization problem in the 
Kac-Moody group recalled in Section 5.1, eq. (pO|) . The point is that given the vacuum wave function ^0 and 
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any element g = g^^g+ of the afRne SL{2, R) Kac-Moody group with the factorisation g± e exp {B± (Ck), 
then the wave function 

* = (*o.9*o^)_*o.9l^ - (*o5*o^)+*o3+^ 

is a sohition of a compatible auxiliary linear system ( [7l| ) . 

Given a wave function ^, the original fields Q±,P± and (p., a, solutions of the Einstein equations, are 
then reconstructed by evaluating matrix elements of the wave function. This is done with the help of vertex 
operators. Choosing as in ref.[[l6[ highest weight representations of the Kac-Moody group in which the group 



elements g can be written as products of vertex operators one recovers formulas (eq.(Mp5l36)) for the Lax 
connection, and more interestingly formulas for the elements of the metric (eq.(plr). (More details may be 
found in ref . pq| . ) We mark that in formula ( pl| ) all the coordinates' dependence is contained in the vacuum 
wave function "^q. 

Note that in this approach the Lax connection does contain any space-time dependance spectral param- 
eter. There is no moving spectral parameter. In particular the poles, which are the w-argument of the 
vertex operators, are fixed in contrast with other methods based on the Belinskii-Zakharov approach. The 
space-time dependence comes back when we conjugate the vertex operators with ^o, cf eg. (p7|) . This is 



manifest in the definition (70) of the moving pole which may be rewritten as: 



So the Virasoro algebra appears as a way to encode the coordinates' dependence of the moving poles. 

8 Appendix B: Vertex operator expectation values. 

Here we have gathered a few formulas for the vertex operator expectation values. These are computed using 
Wick's theorem: 

_ / \ Up-Uq/2 

To evaluate the tau-function (B3) or its dual one needs to compute expectation values such as 

■m n 

(n^«.('^p)-n(i+*^^-^2(A^,))> 

This is done using Wick's theorem (|72|). For example 

n n / \ 2 

The determinant formula for the tau-functions then follows from the Cauchy determinant formula 

det f^^) = TT ("^^^^^ ' 



Uj Hj ^ 



To evaluate the tau-functions ( p5| , [36D one needs to know 






i<j ^ ' ' -^ ^ j 

/ \ Ui-Uj /2 



from which equations (p9|) follow 
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9 Appendix C: Some details about the examples 

Here we have gathered a few formulas and identities that can be helpful in the computations used in the 
examples. First, let introduce the two variables ^ and w by 

C - I (V(l - 2 + P)(l -z-p) + ^il + z + p){l + z-p)^ 



= ^ (V(l - 2 + P)(l -Z-p)-y/il + Z + p){l + Z-p) 



Reciprocally 



z — —w^ 

The domain of definition of p and z implies that \w\ < £, < 1. 

We use the notation Hz and Xz for the variables /i and X associated to z as defined in eq.(0). The 
singular values of z we have chosen, eg. z — ±1 or z = oo, lead to simple expressions for the factor using 
fj,z- For example: 

hm Xz [t^i^j^y = li,, V (t^i^j^y^A 



and 

tor z —)■ +CX3 



The variables p,z for z ~ ±1 are related to the variables i^ and w by 



Aii + 1/Vl + M-i/ 1 + C 

/il — 1\ /l+/i_i\ 1 — ui 



^1 + 1/ \1 — //_i/ 1 + w 



^/^i + p.-i J \1 - w' 
The variables Xz for z — ±1 are related to the variables ^ and w by 

■5—!^ ^^— (dp2 - dz^) = const. — ^ , 

Finally we give the parametrisation used to deduce the Kerr solution from the Chandrasekahar-Xanthopoulos 
one. The relation between the parameters (p, q) and the mass and the angular momentum is 



The relations between the coordinates are 



1^ a 



r- M 



2 



^/M' 



RecaU that i?^ = ,,2 ^ ^2 pQg2 6) and D = r^ - 2Mr + a^. 
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